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ABSTRACT 

Taking advantage of the unprecedented combination of sensitivity and angular resolution af- 
forded by the Herschel Space Observatory at far-infrared and submillimeter wavelengths, we 
aim to characterize the physical properties of cold dust within nearby galaxies, as well as the 
associated uncertainties, namely the robustness of the parameters we derive using different 
modified blackbody models. For a pilot subsample of the KINGFISH key program, we per- 
form two-temperature fits of the Spitzer and Herschel photometric data (from 24 to 500 /im), 
with a warm and a cold component, both globally and in each resolution element. We compare 
the results obtained from diff'erent analysis strategies. At global scales, we observe a range of 
values of the modified blackbody fit parameter fic (0.8 to 2.5) and T^ (19.1 to 25. IK). We 
compute maps of our modelling parameters with fic fixed or treated as a free parameter to test 
the robustness of the temperature and dust surface density maps we deduce. When the emis- 
sivity is fixed, we observe steeper temperature gradients as a function of radius than when it is 
allowed to vary. When the emissivity is fitted as a free parameter, barred galaxies tend to have 
uniform fitted emissivities. Gathering the parameters obtained each resolution element in a 
Tc-fic diagram underlines an anti-correlation between the two parameters. It remains difficult 
to assess whether the dominant eff'ect is the physics of dust grains, noise, or mixing along the 
line of sight and in the beam. We finally observe in both cases that the dust column density 
peaks in central regions of galaxies and bar ends (coinciding with molecular gas density en- 
hancements usually found in these locations). We also quantify how the total dust mass varies 
with our assumptions about the emissivity index as well as the influence of the wavelength 
coverage used in the fits. We show that modified blackbody fits using a shallow emissivity (J3 
< 2.0) lead to significantly lower dust masses compared to the 13 < 2.0 case, with dust masses 
lower by up to 50 % if =1-5 for instance. The working resolution aff'ects our total dust mass 
estimates: masses increase from global fits to spatially-resolved fits. 
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1 INTRODUCTION 

The interstellar medium (ISM) plays a key role in the star forma- 
tion history of galaxies. While stars evolve and die, they re-inject 
dust and gas in the ISM that will then be processed to form a new 
generation of stars. Dust is formed in the envelopes of late-evolved 
stars and is transformed in supernova shock waves or stellar winds 
by violent processes such as destruction, vaporization, sputtering or 
erosion ( Jones et al.|T994l[T996l|Serra Diaz-Cano & Jones|2008| ). 
The size of dust grains also increases through processes like ac- 
cretion or coa gulation accretion (|Stepnik et al.||2003[ Dominik & 



ITielensI 19971 IDominik & Dullemond|2008||Hirashita & Omukai 
|200 9l). They participate in the gas thermal balance and screen parts 
of the ISM through the absorption of starlight and subsequent emis- 
sion of infrared (IR) photons. 

Dust grains contribute to the chemistry of the ISM. Molec- 
ular hydrogen predominantly forms on grains, and molecules can 
freeze out on grain surfaces and undergo additional transformation 
( |Hasegawa & Herbst 1993l|VidaH et al.|2004||Cazaux et al. 2005 ). 
Atomic and molecular gas are commonly traced by Hi and CO ob- 
servations. CO can be a poor tracer of molecular gas, especially in 
low-metallicity environments, due to the dependency of the CO- 
to-H2 conversion factor with density and metallicity ( Wilson 1995 , 
Taylor et al.|[T998l [Israel 2000, Leroy et al. 2005, Wolfire et al.. 



2010| ). Weknow that dust and gas are closely tied in the ISM. This 
suggests that used with Hi and CO, dust could also be used as an 
additional tracer of the gas. 

It is crucial to yield an exhaustive inventory on the range and 
distribution of dust temperatures as well as the radiative heating 
sources of dust grains. IRAS data combined with ground-based ob- 
servations (SCUBA on JCMT for instance, observing at 450 and 
850 yL/m), have first led to the detection of the cool phases of dust 
in nearby galaxies ( [Dunne et al.|2000[ [Dunne & Bales |200T| [James | 
[et al.|2002| . The Infrared Space Observatory (ISO) and the Spitzer 
Space Telescope (with a better sensitivity) then enabled us to re- 
solve in detail the warm and cool dust reservoirs and more accu- 
rately model the Spectral Energy Distributions (SED) of nearby 
objects up to 160 fim ( jDraine et al.|2007j ). 

Several issues still have to be investigated. Numerous stud- 
ies thus highlighted the necessity of observing galaxies longward 
of 160 fim at high resolution to properly sample the submillime- 
ter (submm) regime of these SEDs and better characterise the total 
temperature range and dust masses (e.g. [Gordon et al.|2010| ). More- 
over, IR emission is known to be a good tracer of the embedded star 
formation in galaxies ( Kennicutt|1998|[Calzetti|2007|[Calzetti et al.[ 
|2010| ). Nevertheless, the amount of IR emission not directly linked 
with the ongoing star formation is still poorly quantified. To ad- 
dress this issue, we need to understand the heating sources of the 
diff'erent grain populations: dust heating by the young stellar pop- 
ulations in Hii and photodissociation regions (PDR), UV photons 
escaping from those regions or old stellar population. 

The Herschel Space Observatory ( jPilbratt et al.|2010] ) is cur- 
rently observing the FIR/submm wavelength range to probe the 
coldest phases of dust that were poorly unexplored. Its onboard 
photometers R\CS (from 70 to 160 yt/m) and SPIRE (from 250 to 
500 jim) enable the imaging of nearby galaxies at a resolution and 
sensitivity never reached to date by a space telescope. The Planck 
satellite, launched by the same rocket, performs an all- sky survey at 
nine wavelengths between 350 jmi and 1 cm, complementing Her- 
schel coverage, at a lower resolution. [Planck Collaboration|p01 la| 
studies have suggested that some nearby galaxies are both luminous 
and cool (presence of dust with T<20K), with properties compa- 



rable to those of galaxies observed with Herschel at higher red- 
shifts ( [Rowan-Robinson et all2010 , MagneUi et al.|2012| ). Using 
Herschel observations, [Engelbracht et al.| ( 2010) compute separate 
SEDs for the center and the disks of nearby galaxies and find that 
the cool dust is on average 15% hotter in the center than in the disk 
and correlate this eff'ect with morphological type and bar strength. 
Using a radial approach, [Pohlen et al.[ ( [2010) show that SPIRE sur- 
face brightness ratios - used as a proxy for cold dust temperatures 
- seem to decrease with radius in spirals, and conclude that dust 
in the outer regions is colder than in the centre of the galaxies. 
[Skibba et al.| ( |201l| ) also investigate the relation between dust heat- 
ing sources and the star formation activity with metallicity and the 
morphological types of their galaxies. On a more resolved scale, 
,Bendo et al., ( ^2012j derive PACS and SPIRE ratio maps of spi- 
ral galaxies. They find that those ratios seem to be preferentially 
correlated with the global starlight, including evolved stars, than 
with emission from star forming regions, results consistent with 
those obtained by Boquien et al. (2011 ) for the spiral local group 
galaxy M33 or by Auld et al., ( ,2012 ) for extragalactic dust clouds 
of Centaurus A, among others. This supports the idea of a dif- 
ferent source of heating for the cold (15-30 K) dust than that of 
the warm dust. Using Herschel observations of the dwarf irregular 
galaxy NGC 6822, [Galametz et aL] ( [2010| ) show on the contrary that 
SPIRE surface brightness ratios are strongly correlated with the 24 
fim surface brightness, a commonly used tracer of the star forma- 
tion activity. This might indicate that star forming regions are non 
negligible heating sources of cold dust in that object and suggest 
that the source of the dust heating probably diff'er with the environ- 
ment and vary from normal spirals to low metallicity starbursts. 

Observations of galaxies with Herschel and from the ground, 
have also suggested that the properties of cold dust grains could dif- 
fer from those often used in SED models such as blackbody modi- 
fied with a beta-2 emissivity law or use of graphite grains to model 
carbon dust. In metal-poor objects in particular, using graphite 
grains can lead to unphysical gas-to-dust mass ratios (G/D) com- 
pared to those expected from the metallicity of the galaxy ( Me ixner| 
|et al. [20 10) [Galametz et al.[20T0| . Moreover, observations of dwarf 
galaxies in submm often lead to the detection of excess emission 
compared to extrapolations from Spitzer-hdiS^d fits ( Dumke et al.J 
2004, GaUiano et al.'2003 2005 , Bendo et al. 2006, Marleau et al.] 
2006 , Galametz et al. 2009 , O'Halloran et al. 2010 , Bot et al. 20f0]). 
The origin of this excess is still an open issue and several explana- 
tions have been investigated so far: variations of the emissivity of 
the dust grains with temperature, "spinning dust" emission, cold 
dust etc. Many investigations still need to be carried out to disen- 
tangle between those diff'erent hypotheses. 

We have obtained photometric Herschel/FACS (70 to 160 fim) 
and Herschel/SFIRE (250 to 500 yum) observations of nearby galax- 
ies as part of the open time key programme KINGFISH (Key In- 
sights on Nearby Galaxies: A Far-Infrared Survey with Herschel, 
jKennicutt et al.|20l"T] ). The KINGFISH sample comprises 61 galax- 
ies probing a wide range of galaxy types (from elliptical to irregular 
galaxies), metallicity (7.3 < 12+log(0/H) < 9.3) and star-forming 
activity (lO'^ < SFR < 7 yr-^. The KINGFISH galaxies are 
located between 3 and 31 Mpc, leading to ISM resolution elements 
of 0.5 to 5.4 kpc at SPIRE 500 jum resolution {36''). We refer to 
jKennicutt et al.[ ( [20TT] ) for a detailed description of the sample and 
the observation strategy. 

We aim to benefit from the good coverage and resolution of 
Herschel to sample the complete thermal dust emission spectrum 
and probe the spatially-resolved properties of the cold dust grain 
population. We refer to the work of [Aniano et al.[ ( [20T2| (hereafter 



Mapping the cold dust properties with Herschel 3 



Table 1. Galaxy Data 



Galaxy 


Optical 




5q 


Major diam. 


Minor diam. 


Distance 


12+log(0/H) 


SFR 






Morphology 


(J2000) 


(J2000) 


(arcmin) 


(arcmin) 


(Mpc) 






(Moyr-i) 


(X IQP Lq) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




(8) 


(9) 


(10) 


NGC 337 


SBd 


OOh 59m 50.7s 


- 07d 34' 44" 


2.9 


1.8 


19.3 


8.18 


8.84 


1.30 


12.0 


NGC 628 


SAc 


Olh 36m 41.8s 


15d47' 17" 


10.5 


9.5 


7.2 


8.35 


9.02 


0.68 


8.0 


NGC 1097 


SBb 


02h 46m 18.0s 


- 30d 16' 42" 


9.3 


6.3 


14.2 


8.47 


9.09 


4.17 


45.0 


NGC 1291 


SBO/a 


03h 17m 19.1s 


- 41d 06' 32" 


9.8 


8.1 


10.4 


8.52 


9.20 


0.35 


3.5 


NGC 1316 


SABO 


03h 22m 41.2s 


- 37d 12' 10" 


12.0 


8.5 


21.0 


8.77 


9.52 




8.0 


NGC 1512 


SBab 


04h 03m 55.0s 


- 43d 20' 44" 


8.9 


5.6 


11.6 


8.56 


9.11 


0.36 


3.8 


NGC 3351 


SBb 


lOh 43m 57.5s 


lid 42' 19" 


7.4 


5.0 


9.3 


8.60 


9.19 


0.58 


8.1 


NGC 3621 


SAd 


llh 18m 18.3s 


- 32d 48' 55" 


12.3 


7.1 


6.5 


8.27 


8.80 


0.51 


7.9 


NGC 3627 


SABb 


llh20m 13.4s 


12d 59' 27" 


9.1 


4.2 


9.4 


8.34 


8.99 


1.7 


28.0 


NGC 4826 


SAab 


12h56m 42.8s 


21d40' 50" 


10.0 


5.4 


5.3 


8.54 


9.20 


0.26 


4.2 


NGC 7793 


SAd 


23h 57m 50.4s 


- 32d 35' 30" 


9.3 


6.3 


3.9 


8.31 


8.88 


0.26 


2.3 



Galaxy Name 
Morphological type 
J2000.0 Right ascension 
^"^^ J2000.0 Declination 

Approximate major diameter in arcminutes 
Approximate minor diameter in arcminutes 
Distance in megaparsecs 

Mean disk oxygen abundance from jMoustakas et ar]|2010) 
Star formation rate in solar masses per year from Table 1 o f |Calzetti et all ( |2010t 
Total infrared luminosities in the 3-to-l 100 fim in units of 10^ solar from |Kennicutt et ar]j201 1) 



[A12a]) and Aniano et al. 2012b (in prep, full KINGFISH sam- 
ple, hereafter [A 12b]) for an analysis of other dust properties (PAH 
contribution to the total dust mass, distribution of the starlight in- 
tensity, dust mass surface density, among others) in the KINGFISH 
sample using the Draine & Li (2007 ) dust models. The two stud- 
ies are complementary. The following paper particularly focuses 
on investigating the distribution of cold dust temperatures and dust 
mass surface densities in our objects. The SED model applied in 
this work (two modified-blackbody fit) allows for variations in the 
emissivity index of the dust grains and enables us to test the robust- 
ness of applying diff'erent assumptions of emissivity index to our 
extragalactic objects and quantify how this aff'ects the temperatures 
and masses derived. 



This paper is structured as follows. In §2, we describe the Her- 
schel observations and the convolution techniques we apply to de- 
grade images to the lowest resolution (that of SPIRE 500 fim). In 
§3, we discuss the integrated properties of the cold dust (temper- 
ature and emissivity). We also study how uncertainties in the flux 
measurements can lead to degeneracies between parameters. In §4, 
we produce maps of the dust properties of our galaxies and discuss 
the choice of a fixed or free emissivity parameter on the tempera- 
ture and emissivity maps obtained. We also study the distribution 
of the dust mass surface densities in our objects. We finally analyze 
how the total dust mass depends on the type of modified blackbody 
model we apply, the wavelength coverage and the working resolu- 
tion. 



2 OBSERVATIONS AND DATA PROCESSING 
2.1 The sample 

The data taken with the Herschel Space Telescope now enable us 
to properly constrain the cold dust properties in nearby galaxies. 
In this paper, we focus our study on 1 1 galaxies, 9 spirals and 2 
lenticulars (NGC 1291 and NGC 1316) observed as part of the 
Herschel open time key programme KINGFISH. Those galaxies 
were primarily selected because they have also been mapped at 870 
yL/m with the LAB OCA instrument located on APEX in Chile, data 
that will be analysed in a second paper in this series. Our galax- 
ies 1) are observable from APEX (6o <+21°) 2j are large enough 
to use the LAB OCA mapping mode (the field of view of the com- 
plete array is 11^4) 3) have diff'use emission bright enough to be 
detected with LAB OCA. Our sample is thus biased toward rela- 
tively high surface brightness normal spirals. Albrecht et al. (in 
prep) describe the LAB OCA observations and data reduction. We 
will combine LAB OCA data with the results obtained in this pa- 
per to investigate possible 870 yum submm excess in our objects 
(Galametz et al., in prep). The properties of our sample are sum- 
marized in Table [T] Morphologies and sizes as well as references 
and methods used to derive distances, and the total infrared lumi- 
nosity (Ltir) are detailed in |Kennicutt et al.| ( |2011| ). The SFRs are 
derived using a combination of Ka and 24 /urn integrated luminosi- 
ties and taken from |Calzetti et al.l ( |2010| ). This formula is used in 
§4. 1 of this paper. Two values are given for the oxygen abundance, 



one coming from the empirical calibration of Pilyugin & Thuan 
( 2005 ) (PT), the other from the calibration of Kobulnicky & Kew 
|ley| ( |2004j ) (KK) as discussed in details in Moustakas et al.ipoio] ). 
The sample covers diff'erent morphologies and a wide range of bar 
strengths. It contains isolated targets as well as galaxies with dis- 
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Figure 1. Continued. 



galaxy NGC 1097, was observed and studied during the Science 
Demonstration of Herschel. It hosts a prominent starbur sting ring 
with a radius of ~ 900 pcJHummel et al. 1987 ). We refer to |Sand-| 
[strom et al.| ( [20T0l l and |Beirao et al. ( 2010) for PACS photomet- 
ric and spectroscopic studies of this circumnuclear ring. We also 
refer to Hinz et al. (submitted to ApJ) for a radial study of tem- 
peratures in NGC 1291, another galaxy of this sample. We benefit 
from the good spatial resolution of Herschel to resolve with SPIRE 
500 jim (FWHM = 36 arcsec) structures of ~700 pc for the closest 
galaxy of our sample NGC 7793 and ~4 kpc for the furthest galaxy 
NGC 1316. Three-color images (24 / 100 / 250 iim) of the galaxies 
are provided in Fig.[T] 

2.2 The Herschel observations 

The KINGFISH programme observed their 61 nearby targets with 
the PACS and SPIRE instruments between December 2009 and 
June 2011. A description of the data processing steps are given in 
[Engelbracht et al.| ( [20T0|j Sandstrom et al.|([20T0l) and the KING- 
FISH overview paper of pennicutt etalj(^2011j). A few details are 
summarized in the next sections. 



2.2.7 PACS observations 

The PACS photometer observes at 70, 100 and 160 fim with 
FWHMs of 572, 777 and 12'' respectively ( |Poglitsch et aL|2010| ). 
Observations of KINGFISH galaxies with this instrument were 
obtained with 15' -long cross-scans (perpendicular scans) with a 
medium scan speed of 20" s~^ PACS calibration uncertainties are 
~5 % according to Version 4 of the PACS Observer's Manuaj^ The 
data reduction is performed from raw data to Level 1 data with 
the Version 5 of the Herschel Interactive Processing Environment 
(HIPE) using the standard procedure with additional second-level 
deglitching steps. We use the Scanamorphos technique to process 
the data from Level 1. Scanamorphos is an IDL software that sub- 
tracts the brightness drifts caused by the low-frequency noise us- 
ing the redundancy built in the observations. We refer to |Rous-| 
|sel | (|2012|) for a complete description of the software. PACS data 
emerge from the pipeline with a pixel size of 174, V!l and 2V85 
at 70, 100 and 160 jdm respectively. We applied the needed correc- 
tion factors to obtain the latest version (v6) of the responsivity cal- 
ibration. The root-mean- square (rms) noise levels are [0.14-0.18], 
[0.24-0.35] and [0.74-1.9] mJy arcsec-^ at 70, 100 and 160 fim re- 
spectively for our subsample. The Scanamorphos pipeline aims to 
preserve low surface brightness diff'use emission. Nevertheless, the 
PACS instrument is less sensitive than the MIPS instrument to low 
surface brightness emission. Our analysis of the resolved proper- 
ties of our objects will be restricted to pixels for which we have a 
3-sigma detection in all the bands. 



turbed outskirts resulting from past interactions such as NGC 3627. 
It also contains mergers such as NGC 1316 (Fornax A), an early- 
type radio galaxy that could be interacting with NGC 1317, also 
detected in our Herschel field of view. NGC 1316 possesses a core- 
jet lobe structure ( |Ekers etaL] 1983|) and a low-luminosity X-ray 
active galactic nucleus ( Kim & Fabbiano||2003| ). We discuss how 
this could aff'ect our results later in the paper. The star formation 
rates (SFR) of our sample range from 0.3 to 8 M© yr~^ and the 
l^TiR covers an order of magnitude, from 2.3 x 10^ to 4.5 x 10^^ 
Lq. Our sample thus spans a wide range of environments and phys- 
ical conditions to which dust is exposed. One of them, the barred 



2.2.2 SPIRE observations 

The SPIRE instrument produces maps at 250, 350 and 500 //m, 
with FWHMs of 18'', 25'' and 36" respectively ( [Griffin et al.|20T0l ). 
Observations of KINGFISH galaxies with this instrument were ob- 
tained in scan mode. The data reduction was performed from raw 
data with the Version 5 of HIPE. Details on the data reduction are 
given in the KINGFISH Data Products Delivery User's Guide 

^ http://herschel.esac.esa.int/Docs/PACS/html/pacs_om.html 
^ http://herschel.esac.esa.int/UserReducedData.shtml 
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Table 2. Global IRAS, MIPS, PACS and SPIRE flux densities (in units of Jy) 
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54.5±10.9 


115±23.0 


11.5±2.3 


9.2±1.8 


12.6±2.5 


41.1±8.2 


82.1 + 16.4 


137±27.3 


81.7±16.3 


54.1±10.7 


PACSIOO^ 


19.5+1.0 


74.0+3.7 


116±6.0 


12.8+0.7 


9.3±0.5 


13.8+0.7 


46.1+2.3 


94.4+4.7 


179+9.0 


95.7±4.8 


65.8+3.3 


MIPSI6O2 


20.1±2.4 


112±17 


154±19 


26.3+4.1 


12.6±1.8 


19.6±3.1 


56.9±8.9 


139±17 


216±28 


85.8±13 


107±17 


PACS1604 


19.6+1.0 


116±6.0 


134±7.0 


20.3+1.1 


11.5+0.6 


18.7+1.0 


55.1+2.8 


128±6.0 


202+10 


94.1+4.7 


91.2+4.6 


SPIRE250'^ 


9.8±0.7 


65.5±4.7 


72.2±5.1 


15.9±1.1 


4.8±0.4 


15.6±1.1 


32.4±2.3 


71.2±5.1 


96.7±6.9 


42.4±3.0 


56.3±4.0 


SPIRESSO"^ 


4.4±0.3 


30.6±2.2 


30.8±2.2 


8.0±0.6 


2.1±0.2 


8.7±0.6 


13.7±1.0 


31.7±2.3 


37.6±2.7 


16.4±1.1 


28.4±2.0 


SPIRESOO"^ 


1.9±0.1 


13.3±1.0 


12.6±0.9 


3.5±0.3 


0.8±0.1 


4.2+0.3 


5.3+0.4 


14.6+1.0 


14.4+1.0 


6.3+0.4 


13.9±1.0 



Parameters of the photometric apertures used to derive Herschel flux densities. 2a and 2b : lengths of the major and minor axes; a: position angle of the 
apertures major axis measured east of north. 
^ from |Dale et alH2007i . 

^ from the SCANPI tool and the HIRES atlas (see text). 
4 from |Dale et a'l1 ( |2012} . 



The SPIRE maps were built with a nearest-neighbor projection on 
sky and averaging of the time ordered data. SPIRE maps emerge 
from the pipeline with a pixel size of 6" , 10'' and 14'' at 250, 350 
and 500 jim respectively We use the SPIRE beam areas quoted in 
the SPIRE Observer's manual, namely 423, 751 and 1587 square 
arcseconds at 250, 350 and 500 fim, with ~1% uncertainty The rms 
noise levels are [0.5-0.75], [0.24-0.34] and [0.14-0.21] MJy sr-^ at 
250, 350 and 500 jim respectively for our subsample. Calibration 
uncertainties are estimated to be ~7% for the three wave bandfl 



2.3 Background galaxies and sky subtraction 

The background galaxies subtraction and sky subtraction steps are 
similar to those applied in ' Dale et aE] ( [2012 ') to derive the global 
photometry for the KINGFISH sample. We identify the bright 
background galaxies that could contaminate the photometry and 
"remove" them by replacing the data with random values with sim- 
ilar noise characteristics. The removal of these background objects 
does not significantly aff'ect the total flux estimates (<1%). The 
companion of NGC 1316 (NGC 1317) is also removed from the 
Herschel fields of view to avoid strong contamination in the chosen 
photometric aperture. A global sky subtraction is already applied 
on the PACS/SPIRE maps during the reduction process. However, 
this value was estimated on the entire map. We determine a refined 
local sky using apertures randomly distributed in the background. 
We choose the total "sky" area to be greater than the area we choose 
to perform the photometry (see Table[2|. We do not expect Galactic 
cirrus to contaminate any of the galaxies of our sample. Indeed, the 
lowest Galactic latitude for the sample is +26° for NGC 3621. 



http://herschel.esac.esa.int/Docs/SPIRE/html/spire_om.html 



2.4 Convolution to lower resolution maps 

To compare the fluxes at diff'erent wavelengths, we need to con- 
volve every observation to a common resolution. We refer to Ani-| 
,ano et al. ( 2011) for details on how to construct convolution kernels 
between Spitzer and Herschel bands. The main steps are similar to 
the technique described in [Gordon et al.| ( [2008] ). We review the ba- 
sic steps here. 

In-flight Point Spread Functions (PSFs) are available for Her- 
schel. An ideal convolution kernel K(x,y) which transforms an orig- 
inal PSF (PSFJ to a lower resolution final PSF (PSFj) is defined 
as: 



K{x,y) = FT- 



FT[PSFf(x,y)] 



FT[PSFo(x,y)] 



(1) 



To attenuate the high-frequency noise in the original PSF, a 
filter is first applied to the FT[PSF^(x,y)]: 



Filter{(jS) - 



1 



^ 1 1 + cos 



for 0) < ojiow 
for ojhigh < 



where cohigh is the cutoff" spatial frequency and coiow- 0.7 cohigh- (^Mgh 
is appropriately chosen to reach a high accuracy but minimizes the 
noise linked with the process. The final kernels are circularized 
to produce rotationally symmetric kernels, thus independent of the 
relative orientation of the spacecraft. 

We convolve MIPS, PACS and SPIRE images to the SPIRE 
500 jim images using the previous kernels and the IDE procedure 
convJmage. We finally project the maps to a common sampling 
grid (that of the SPIRE 500 jjm maps) with the IDE Hastrom ^xo- 
cedure. [Aniano et aL] { |201 1 ) caution about convolutions from MIPS 
160 to SPIRE 500 iim. The deconvolution of an image to higher 
resolution can indeed be risky if the image contains noise and arte- 
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facts. We compare the resolved PACS 160 jdm flux densities (pixels 
about a 3 -cr in the band) for the whole sample with the resolved 
MIPS 160 jim flux densities (maps convolved to SPIRE 500 res- 
olution). We find a very good agreement between the two values 
and therefore decide to use both in the modelling. As mentioned 
in [A 12a], among others, slight discrepancies can be observed be- 
tween MIPS and PACS fluxes. We know for instance that the PACS 
instrument is less sensitive to diff'use emission than MIPS. Using 
both MIPS and PACS (at 70 and 160 iim) data enables us to be 
more careful and conservative in our modelling. The model will 
indeed try to find a compromise between the two values. 



3 INTEGRATED PROPERTIES 
3.1 Naming convention 

The SED peak emission of nearby galaxies spans a relatively nar- 
row range of dust temperatures, from ~15 to 40 K (|Bendo et al. 



NGC 0337 



2003 , Stevens et al .|2005[|DaIe et al.|2005[|Bosem et aL|2010[|Dale 
et al. 2012 for instance). In this study, we decide to use an isother- 
mal fit to describe the cold dust component and to call our param- 
eter Tc the "cold dust temperature". It is nevertheless important to 
keep in mind that the temperature distribution of the cold dust phase 
is more complex than a delta function, with a range of cold grain 
temperatures along the line of sight. The "cold dust temperatures" 
derived in this paper (and in other studies using similar fitting pro- 
cedures) should be considered in reality as an approximation of the 
cold temperature distribution in our galaxies. Moreover, an emis- 
sivity index is, by definition, the power-law index describing the 
variation of dust opacity with wavelength, so an intrinsic property 
of the grains. Calling our parameter pc the "cold dust emissivity 
index" is a commonly-used expression. Our parameter pc is more 
accurately the "eff'ective emissivity index", namely a power-law in- 
dex resulting from the sum of diff'erent grain populations, each of 
them having their own intrinsic emissivity index. The eff'ective and 
intrinsic emissivity indices are identical in the isothermal case. 

3.2 The model 

We first want to derive integrated dust properties of our sample to 
investigate the variations of the cold dust temperature and emis- 
sivity between objects. A combination of data from 24 to 500 jim 
provides a complete sampling of the warm and cold dust thermal 
emission of our galaxies from which global dust properties can be 
deduced. The submm data are, in particular, crucial to constrain the 
coldest phases of dust that could account for a large amount of the 
dust mass. They also enlighten us as to how the properties of cold 
dust vary from galaxy to galaxy. 

We use the gl obal PACS and SP IRE fluxes of KINGFISH 
galaxies derived by |Dale et aT] p012| ). Their elliptical apertures 
(Table |2]) are similar to those used to perform the Spitzer photom- 
etry (for SINGS in particular) and encompass the total emission 
of the galaxy at each wavelength. Global fluxes are measured 
on original resolution images. We complement the global SED 
coverage with data from Spitzer/MIFS and IRAS. We use the 
global flux densities of our galaxies for Spitzer/MIFS bands 
derived by [Dale et al.| ( [2007} . We obtain the IRAS global flux 
densities at 25, 60 and 100 jum (accurate to ~20%, when available) 
through the SCANPj^tool (IRAS Scan Processing and Integration 





0.1 Ir 




|l . -I-Y--'---A--^->---0-| 



-0.1 I 
-0.2 I 



100 

Wavelength (/xm) 




100 
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100 1000 
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Figure 2. Global SED of our galaxies performed using a two-temperature 
fit. The emissivity index fi^ of the warm component (dashed dotted line) 
is fixed at 2.0 and the emissivity index j3c of the cold component (long 
dashed line) is a free parameter. IRAS (asterisks), MIPS (squares), PACS 
(downward triangles), SPIRE (upward triangles) data are overlaid. SCUBA 
850 yum flux densities are indicated with empty circles but not used in the 
fit. If the error bars are not shown, errors are smaller than the symbols. The 
fitted /3c, and temperatures are given in Table [3] The bottom panel of each 
plot indicate the residuals from the fit. See §3.1 for details on the modelling 
technique. 



^ http://scanpiops.ipac. caltech.edu:9000/appHcations/Scanpi/index.html 
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Figure 2. Continued. 
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Table 3. Best fit parameters 







He — ^ 






R free 








WO/zm, warm/ ^lOfxm, tot 






-LvO/im, warm/ ^lOfim, tot 


^w 


He 




NGC 337 


0.58 


51.4 ± 1.7 


21.3 ±0.8 


0.35 


55.7 ± 14.3 


1.6 ±0.2 


25.2 ± 2.4 


NGC 628 


0.53 


55.9+ 1.8 


19.0 ± 0.4 


0.52 


56.0 ± 5.3 


2.0 ± 0.2 


19.1 ± 1.3 


NGC 1097 


0.57 


55.8 ± 1.6 


20.7 + 0.6 


0.57 


55.9 ± 2.4 


2.0 ± 0.2 


20.7 ± 1.4 


NGC 1291 


0.57 


56.2 ± 2.0 


17.6 + 0.4 


0.36 


60.8 ± 10.2 


1.5 ±0.2 


20.4 ± 1.4 


NGC 1316 


0.32 


58.9 ±8.1 


22.4 ± 0.8 


0.52 


55.2 ± 11.2 


2.5 ± 0.4 


19.4 ± 2.4 


NGC 1512 


0.81 


49.9+ 1.1 


16.8 ± 0.5 


0.38 


56.4 ± 16.8 


0.8 ± 0.2 


25.1 ±2.9 


NGC 3351 


0.51 


57.6 + 2.2 


20.4 ± 0.5 


0.53 


57.4 ± 3.3 


2.1 ±0.2 


20.0 ± 1.2 


NGC 3621 


0.54 


54.9+ 1.6 


19.7 ± 0.5 


0.47 


55.9 ± 5.3 


1.8 ±0.2 


20.7 ± 1.4 


NGC 3627 


0.32 


59.2 ±6.1 


22.4 ± 0.6 


0.47 


55.8 ± 5.6 


2.3 ± 0.2 


20.2 ± 1.4 


NGC 4826 


0.21 


59.1 ± 13.5 


23.7 ± 0.7 


0.24 


57.0 ± 16.4 


2.1 ±0.2 


23.1 ± 1.6 


NGC 7793 


0.59 


51.9 ± 1.7 


18.7 ±0.5 


0.17 


63.4 ± 19.5 


1.2 ±0.2 


24.4 ± 1.9 



Note : The temperatures and emissivity indexes given in this table are the median values of the Monte-Carlo distributions (see §3.2). 



tool) and the HIRES (IRAS High Resolution Image Restoration, 
ISurace et al.p004| atlas. We provide the MIPS, PACS and SPIRE 
flux densities in Table [2] MIPS and PACS fluxes are in very 
good agreement. We refer to [Dale et aT] ( |2012| for quantitative 
comparisons between MIPS and PACS global fluxes at 70 and 160 
fim. 

Cold dust temperatures are often derived using single modified 
blackbody (MBB) fits. Nevertheless, a galaxy contains a range of 
dust temperatures that a single MBB will not be able to capture. The 
warm dust component, in particular, spans a narrow range of dust 
temperatures, from ~40 to 100 K and contributes to a non negligi- 
ble fraction of the 70 fim emission. As a result, using a single MBB 
model fit usually leads to an overestimation of the cold dust temper- 
atures. We note that this effect particularly affects objects showing 
an enhanced 24-70 fim emission linked with stochastically heated 
grain emission, like the Magellanic Cloud for instance fB ot et al^ 
[2004 1 [Bernard et aL||2008|). We thus carry out a two-temperature 
modeling of our galaxies in order to separate the warm and the 
cold component, using the MIPS, PACS and SPIRE flux densities 
Usted in Table[2]along with MIPS 24 and IRAS 25, 60 and 100 
yL/m flux densities. We fit the data dy with a model of the form: 

dy = A, [Ar^By{T^)) + A, [x-^^By{TS) , (3) 

with By the Planck function and the 5 free parameters of our 
model: Pc is called the emissivity index of the cold component, 
Tw and Tc the temperature of the warm and cold component 
respectively and A^ and Ac the overall amplitudes of each dust 
component (see discussion on the naming convention in the next 
section). This equation assumes optically thin emission over our 
wavelength range. As can be observed in Eq[3] we choose to 
fix the emissivity index of the warm dust component in order 
to limit the number of free parameters in our model. We decide 
to fix Pw to 2.0. This is a good approximation of the opacity 
of the graphite/silicate dust models of |Li & Draine| p001| ). We 
note that our choice of p^ does not influence the cold dust 
temperatures. Fixing p^ to 1.5 leads to a small decrease of the 
cold dust temperatures of less than 0.6% if /?c=2.0, and less than 
1.6% if Pc is free. [Tabatabaei et al.| ( |2011) (see also Tabatabaei 
et al., in prep) have synthesized the two MBB approach for 
M33, with Pc free and p^ fixed to 2.0, 1.5 and 1.0 and found that 



Pw-2.0 best reproduces the observed fluxes, reinforcing our choice. 

We compute the expected photometry for the fitted model by 
convolving this model with the instrumental spectral responses of 
the different cameras. We refer to the individual users manuals of 
each instrument for the usual conventions used for each instrument. 
We use the IDL function MPCURVEFIT (Markwardt 2006) to 
perform a Levenberg-Marquardt least- squares fit of the data. We 
weight the data during the minimization in order to take the 
uncertainties on flux measurements into account. 



3.3 Uncertainties on the parameters 

[Shetty et"aL] p009| ) among others show that fits performed solely 
in the Rayleigh- Jeans regime are very sensitive to instrumental 
errors. Fits using the full FIR-to-submm coverage usually lead to a 
significantly better determination of dust properties than fits only 
performed near the SED peak. In order to quantify the uncertainties 
on the parameters driven by the errors of our individual fluxes, we 
perform Monte Carlo simulations for each galaxy. We generate 
500 sets of modified constraints. Each flux randomly varies within 
its error bar, following a normal distribution around its nominal 
value. The absolute calibration is correlated across the three 
SPIRE bands (see the SPIRE 's observers 's manual), namely flux 
densities in the three bands are expected to move in the same 
direction if the absolute calibration would be revised. To account 
for this correlation, we vary the three SPIRE measurements consis- 
tently; i.e., they are not allowed to vary relative to each other. We 
finally apply our two-temperature fit to these 500 modified datasets. 



3.4 Global dust temperature and emissivity index 

Figure |2] shows the global SEDs we obtained. Pc is a free param- 
eter in those models. We overlay the SCUBA data on our SEDs 
when available (indicated with open circles in Figure |2]). Some of 
the SCUBA maps do not fully cover the outer part of the galaxies 
or show a poor detection at the rim of the disks due to a possible 
removal of diffuse emission during the data reduction. Those data 
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Figure 3. Distribution of the cold dust parameters in the temperature-emissivity index plane. 500 fits of randomly modified datasets were performed. The 
green crosses indicate the distribution when the correlation between SPIRE calibration errors is applied, the black crosses when no correlation was taken into 
account. The red dot indicates the median values on T^ and (Sc- Those values are reminded in each panel. The last panel summarizes the trends of all the 
galaxies. 



are not used in the modelling. Bottom panels for each galaxy in- 
dicate the residuals from the fit in logarithm scale. The diflferences 
in shape of the SEDs we obtain highlight the wide range of local 
ISM environments covered by the sample. Discrepancies particu- 
larly appear on the Wien side of the SEDs. We however note that 
the warm component does not significantly afifect the fit to the long 
wavelength cold component. The lOO-to-24 fim flux density ratios 
vary from -17 for NGC 1097 to -37 for NGC 4826, thus from ac- 
tively to quiescent star forming galaxies. Our SEDs mostly peak at 
similar wavelengths, residing between X-\\l yt/m for NGC 337 and 
/l=165 jdm for NGC 1512, in their SEDs. We indicate the contri- 
bution of the warm component to the 70 jim fluxes in Table [3] This 
contribution can be higher than 50% in some objects, legitimizing 
the necessity of two dust components in our simple approach. We 
also note that the relative fluxes at 60 and 70 jjm influences the fits 
as it samples the SEDs where the warm and cold dust overlap. We 
try to fit our data excluding the IRAS 60 yt/m data. This leads to 
a decrease in our cold dust temperature by less than 1.8% if Pc is 
fixed to 2.0 and less than 8.7% if Pc is free. 

We summarize the median values of T^ and pc in Table|3] with 
their errors, namely the standard deviation of the parameter distri- 
bution obtained using the Monte-Carlo simulations. We also indi- 
cate the cold dust temperatures obtained with an emissivity index Pc 
fixed to a standard value of 2.0 for comparison. The median reduced 



ch-squared value is ~1.4 when pc is fixed to 2.0 and ~1.1 when pc 
is free, which means that both models lead to satisfying fits. When 
pc is free, we observe variations of the cold dust emissivity within 
our sample, with values ranging from 0.8 to 2.5 (Table [3]). Those 
values are consistent, within the uncertainties, with the emissivity 
index ranges quoted in the literature ( Hilde brand|1983[pDupac et aT] 
[2003 1 among others). We note that using a free emissivity index 
also increases the uncertainties on the temperature estimates. 

Figure |3] enables us to visualize the distribution of the 500 
runs, results of the Monte-Carlo simulations, for each galaxy in a 
Tc - pc parameter space. In all panels, the green crosses are (Ic^Pc) 
combinations derived from 500 fits taking the correlation between 
SPIRE calibration uncertainties into account. For comparison, we 
also plot with black crosses the (T^, pc) combinations obtained if 
SPIRE data are allowed to vary independently from each other. 
The median values of the temperatures and emissivity indexes are 
annotated in red. We remind the reader that those are the median 
values reported in Table [3] The last panel gathers the individual 
trends. Those trends are obtained by taking the average every 30 
runs. We note that, as expected, taking into account the correla- 
tions between SPIRE calibration uncertainties narrows the ranges 
of values reached by the parameters. Nevertheless, this assumption 
does not change the median values of T^ and pc obtained (less than 
2% diflference). We observe a clear Tc-Pc anti-correlation in all pan- 
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Figure 4. Upper left: MIPS 24 //m image (in MJy sr~i) of NGC 3627. Upper right: Cold dust temperature map of NGC 3627 (in K) derived using /?c=2.0. We 
overlay the MIPS 24 jim contours (smoothed to the resolution of the temperature and beta maps) for comparison. The FWHM of the PSF of SPIRE 500 fim is 
shown by the dashed circle. Lower left: Emissivity index map of NGC 3627 derived with Tc and fic as free parameters. The value ySc=2.0 is indicated on the 
color scale. Lower right: Cold dust temperature map of NGC 3627 (in K) derived with T^ and j3c as free parameters. For each map: North is up, east is left. 
Temperature and emissivity maps for the rest of the sample are available in the appendix of this paper. 



els of Fig. [3] f Shetty et al.| ( |2009| ) caution against the interpretation 
of the inverse Tc-j3c correlation as a physical property of dust in 
the ISM since it could be created by measurement uncertainties. 
Our results clearly confirm that fitting SEDs to noisy fluxes, even 
with data sampling the thermal dust emission from 24 to 500 fim, 
leads to T^-fic anti-correlation. We further discuss this eff'ect in the 
following section. 



4 MAPPING THE COLD DUST TEMPERATURE AND 
EMISSIVITY 

4.1 Dust temperature distribution - /3c fixed 

The spatial resolution of the Herschel instruments now enables us 
to reproduce the previous study on a local basis in order to resolve 
the dust temperature and emissivity within galaxies and directly 
link them with other ISM tracers. We use our two-temperature fit- 
ting approach to derive maps of dust temperature and dust emis- 
sivity in our objects. In order to constrain the warm dust emis- 
sion, we use Spitzer/MIFS images at 24, 70 and 160 yum. These 
maps were obtained as part of the SINGS program (Spitzer Infrared 
Nearby Galaxies Survey; [Kennicutt et al.|2003| ) and reduced using 



the SINGS Fifth Data Delivery pipeline. MIPS, PACS and SPIRE 
maps are then convolved to the SPIRE 500 fim resolution (see §2.4 
for the method) and projected to a common sample grid with a 
pixel size of 18'^ thus half of the FWHM of SPIRE at 500 yum. 
This pixel size corresponds in linear scale to ISM elements rang- 
ing from ~ 340 pc for our closest galaxy NGC 7793 to ~ 1.8 kpc 
for the furthest galaxy of our sample NGC 1316. The emissivity of 
the warm component is fixed to a standard emissivity of 2.0. We 
first use a fixed emissivity index for the cold dust component to ob- 
tain a robust determination of the cold dust temperature distribution 
based on color temperatures. We study the case with a free Pc in the 
following section. 

We exclude pixels with flux densities below a 3-cr detection 
limit in all bandpasses. Uncertainties on the fit parameters are de- 
rived on a local basis using Monte Carlo simulations as described 
in §3.2. We generate 100 sets of modified constraints per pixel and 
run our SED modelling technique on those sets. We average (me- 
dian) the 100 dust temperature values obtained in each pixel and 
derive the final cold dust temperature map. We show the cold tem- 
perature map of NGC 3627 obtained with Pc fixed to a standard 
value of 2.0 in Fig. |4] (upper right panel). We overlay the MIPS 24 
jim contours on the maps to compare the temperature distribution 
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Figure 5. Comparison between the cold dust temperature obtained with our two-MBB models and the cold dust temperature derived from the radiation field 
of the diff'use ISM component (\Jmin) from the |Draine & Li[p007 > models [DL07]. Squares are pixel-by-pixel values (pixel size: 18'0 averaged 4-by-4 for 
clarity. The line indicates a ratio of unity. 



with the star forming regions of the galaxies whose 24 jdm emission 
is a usual good tracer. We show the cold dust temperature maps 
of the whole sample in an appendix in Fig. |A1| (top middle panel 
for each galaxy). We observe a clear radial trend of dust temper- 
atures in our objects. The temperature maxima coincide with the 
center or the star forming regions. The dust grain temperatures de- 
crease in the outer parts of galaxies. This radial trend is consistent 
with results obtained with Herschel observations in nearby galaxies 
(Bendo et al. 2012 , Boquien et al. 201 1 , Engelbracht et al.llOlOl or 
Hinz et al. in prep). 

The source of heating of the cold dust is still poorly un- 
derstood. The emission of old stars, traced by the IRAC 3.6 fim 
observations, smoothly decreases with radius, at our resolution, in 
our objects. A radial trend of dust temperatures may thus suggest 
that the old stellar populations could contribute to the heating of 
the cold dust while the youngest star-forming population is closely 
associated with the warm dust (Calzetti 2007 ). However, the star 
formation activity of our galaxies tends also to decrease radially 
in some of our objects as suggested by the 24 jjm emission, and 
the resolution of our dust temperature maps may not allow us 
to properly disentangle between the different possible heating 
sources. We observe nonetheless that the cold dust temperature 
distribution follows the bars in the galaxies NGC 1097 and 
NGC 3627 or strongly peaks in the nucleus of NGC 1316 where 
the 24 jim emission peaks as well (see Fig. |A1| top left panels), 
suggesting a more significant role of the star forming regions than 
the ISRF in heating the cold dust in those objects. 

[A 12b] perform a pixel-by-pixel SED modeling of the whole 
KINGFISH sample using the |Draine & Lil ( [2007 l) dust models on 
Spitzer + Herschel data up to 500 jim. These models couple a "PDR 



component" reproducing the emission of photodissociation regions 
with a range of radiation fields and a "diflfuse ISM component" 
exposed to a single radiation field intensity U^/^. We note that in 
this model, most of the cold dust is expected to reside in the diflfuse 
component. 

The energy absorbed by a cold grain in this diflfuse phase is: 
QabsWUA^)dAocU^i,, (4) 



with Qabs the absorption efficiency of the grains and U the intensity 
of the radiation field of the diffuse ISM, thus Umin- 

The energy emitted by a cold grain in this diffuse phase is: 

v^BATe,)dvocT'^-P, (5) 



with By the Planck function and Teq the equilibrium temperature in 
the phase, namely the temperature of the cold dust T^. 

Assuming thermal equilibrium conditions, we can write that 
Tabs-^em- Wc Can thus associate the radiation intensity U^/^ and 
the cold dust temperature T^. If we assume that dust grains pos- 
sess an emissivity index /?=2.0 and normalize the temperature to 
the equilibrium dust temperature of the Galaxy (17.5K, |Boulanger| 
fet aL|1996| ), we obtain: 

T,{K) = 17.5 X t/ir^^ = 17.5 X U]^. (6) 

We can thus compare our derived cold temperature maps 
obtained with a fixed emissivity j3c=2.0 with the Umin map obtained 
by [A 12b] using this formula. Dust temperatures have very similar 
distributions in both cases. We plot the pixel-by-pixel comparison 
in Fig [5] We observe that in regions showing a "hotter" cold 
dust temperature (co- spatial with star forming regions or high 
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Figure 6. Dependence of the dust temperatures with the chosen cold emissivity index. Left: Global obtained using Pc = 1. 1-5 and a standard emissivity 
index of 2.0. Galaxies are color-coded. Middle: Median of the pixel-by-pixel obtained for star forming regions. Right: Median of the pixel-by-pixel 
obtained for less active regions. 



starlight intensities), the temperatures obtained in both cases are 
very similar. The points diverge in regions which show a lower 
dust temperature, with colder temperatures deduced using the 
Umin parameter obtained from the [Draine & Li| ( |2007| ) and the 
conversion of Eq. 6. This could be explained by the fact that the 
[Draine &Ti|p007 l model assumes a range of dust temperature and 
thus allows to reach colder dust than a MBB fit with j3c fixed to 2.0. 
Our hypothesis of an emissivity index of to translate U^/„ 

into temperatures slightly influences our results but not enough to 
explain the trend we observe. We refer to [A 12b] for a complete 
description of the Vmin maps for the whole KINGFISH sample. 

We also derive temperature maps with an emissivity index of 
the cold dust component fixed to ySc=1.5 and ySc=1.0. This choice 
does not modify the relative distribution of cold dust temperatures 
in our objects but influences the temperature ranges obtained, scal- 
ing temperatures to higher values than in the = 2.0" case. We 
want to compare this effect in bright star-forming regions as well 
as in more quiescent regions. We separate our pixels into two dif- 
ferent categories based on their star formation rate, using the MIPS 
24 jim band as a qualitative indicator for the SFR and the relation 
from |Calzetti| < [20571 ) : 

SFR(Mq yr-^) = 1.24 X lO-^\L(24iim)(erg s-^)f^^ (7) 

We fix the separation between active and less active regions at 
SFR = 4 X 10~^ Mq yr~^ kpc~^ and plot the variations of dust 
temperatures as a function of the emissivity index in Fig. |6] The 
left panel indicates the dust temperatures obtained from integrated 
fluxes using /3c = 1.0, 1.5 and 2.0. Galaxies are color-coded. The 
errors are derived from our Monte Carlo simulations performed on 
global fluxes (§3.2). The middle panel shows the median values of 
the pixel-by-pixel dust temperatures if we restrict the study to the 
bright star forming regions of the galaxies. The right panel shows 
the median values of the pixel-by-pixel dust temperatures if we re- 



strict the study to the more quiescent regions of the galaxies. For 
the middle and right panels, errors are the standard deviations of 
the distributions of dust temperatures. 

As mentioned earlier, we observe that using a lower emissiv- 
ity index fic shifts the temperature range to higher values. The me- 
dian temperatures of bright star forming regions are systematically 
higher than the median temperatures of more quiescent regions, as 
intuitively expected. With the temperatures increase by up 

to 25% for the bright regions and by up to 21 % for quiescent re- 
gions. With j3c=\.0, the temperatures increase by up to 51% for the 
bright regions and by up to 48 % for less active regions. The choice 
of fic has a major impact on the dust temperatures obtained and 
on the dust masses usually derived from those blackbody fits (see 
§4.3). Our comparison between quiescent and active regions seems 
to indicate that the choice of /3c has a non-linear influence on the 
dust temperatures derived, with a slightly larger impact in bright 
regions. 



4.2 Investigations on emissivity index variations 

The coverage and resolution of Herschel enable us for the first time 
to properly constrain the submillimeter slope of resolved SEDs, 
allowing for more complex analysis of potential variations in the 
dust emissivity at local scales. To investigate these variations, we 
apply the same local MBB fitting technique with both Tc and /3c 
as free parameters. We show the emissivity index map and the 
cold temperature map of the galaxy NGC 3627 obtained with 
those assumptions in Fig. [4] (lower panels). We gather the cold 
dust temperature maps of the whole sample in Appendix Fig. |A1| 
(bottom left and middle middle panel for each galaxy). We overlay 
the MIPS 24 jjm contours as a comparison. The same color scale 
is used for the two cold dust temperature maps. We note that 
the choice of /3c (2.0 or free) strongly influences the cold dust 
temperature ranges obtained for NGC 337 and NGC 7793. We 
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Figure 7. Temperature dependence of the spectral index in our sample. Global estimates are shown with stars. Squares are pixel-by-pixel values (pixel size: 
\W) averaged 4-by-4 for clarity. 



used different color scales for these galaxies. 



Temperatures - Choosing pc as a free parameter modifies 
the temperature distribution in our galaxies. The dust temperature 
distribution is less homogenous in the outer parts of galaxies (in 
NGC 628 or NGC 1097 for instance). This is probably linked 
to the fact that our model becomes very sensitive to noise at 
low- surface brightnesses and was previously observed in |Planck| 
Collabor ation' ('2011b') for the Large Magellanic Cloud and in 
Foyle et al. (2012 ) for M83. When p is used as a free parameter 
in our model, we observe that the distribution of cold dust tem- 
peratures follows a systematic behavior with the morphology of 
the galaxy. Barred spiral galaxies such as NGC 1097, NGC 1291, 
NGC1512, NGC 3351 and NGC 3627 show distributions of their 
dust temperature map similar to those obtained with a fixed Pc, 
i.e. a distribution that follows the star forming regions or that is 
decreasing in the outer part of the galaxies. NGC 3627 presents an 
enhancement of the dust temperature at the end of its bar. Bars alter 
the abundance distribution of the gaseous component in galaxies 
([Friedli et al.|1994| and aff'ect the star formation activity ( |Ellison| 
|et al.||2011| ). Strong concentrations of molecular gas are detected 
at the end of the bar of NGC 3627 ( [Warren et aLl[20T0l ) which 
is probably responsible for the high star formation activity and 
thus for the high dust temperature observed. On the contrary, the 
grand-design spiral NGC 628 or bulge-less objects like NGC 3621 
or NGC 7793 show cold dust temperatures that are homogenous 
throughout the structure of the galaxies. The fact that the cold 
dust temperature distribution does not seem to correlate with any 
dust heating source (neither star forming regions, nor old stellar 
populations or radius) within the galaxies is worrying and ques- 
tions the use of a free pc in the model, at least for this type of galaxy. 



Emissivities - Here again the emissivity index maps seem 
to show a distribution that depends on the morphology of the 
galaxy. Barred spiral galaxies show an homogeneous distribution 
of emissivity index, with median values similar to the pc values 
determined globally. On the contrary, we observe a radial decrease 
of the emissivity index in the non-barred galaxies NGC 628, 
NGC 3621 or NGC 7793 when Pc is free. This radial decrease 
of pc was also observed in M33 by [Tabatabaei et al.| ( |2011| ) or in 
Andromeda by [Smith et al.| ( [20T2] ). We nevertheless note a clear 
inverse relation between temperature and emissivity index in some 
of these objects. We discuss this further in the following section. 
We do not observe a clear trend in NGC 337 and NGC 4826. 
Finally, we obtain low local emissivity values for the galaxies 
NGC 337 and NGC 7793 compared to the standard value of 
2.0 as shown by the emissivity index ranges of each maps. 
This low emissivity indices are consistent with the integrated 
values obtained at global scales. Indeed, we find global emissivity 
indices Pc of 1.6 and 1.2 for NGC 337 and NGC 7793 respectively). 

Tc - Pc anti-correlation - Figure |7] gathers pixels that show 
a detection above our signal-to-noise threshold in each Herschel 
band for our 11 galaxies in a T^ - pc diagram. Global values 
of Table [3] are also overlaid with stars. We estimate the errors 
using the standard deviation of the Monte Carlo simulations of 
uncertainties performed for each pixel. They include calibration 
and sky background uncertainties. An anti-correlation is clearly 
observed between the temperature and the emissivity, with higher 
emissivity values in colder ISM elements. This relation is coherent 
with laboratory measurements (M ennella et al.| |1998[ |Meny| 
|et aL||2007l and references therein) that have reported a strong 
dependence between the temperature and the spectral index of 
the absorption coefficients, with higher p values found for dust 
grains in cold cores. The anti-correlation between T^ and Pc does 
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Figure 8. Dust mass surface density map of NGC 3627 (in M© kpc ^, log scale) derived with /?c = 2.0 (left) and [3c free (right). We overlay the MIPS 24 jim 
for comparison. The dust mass surface density maps for the rest of the sample are available in the appendix of this paper. 



not disappear if we restrict ourselves to the pixels showing the 
smallest uncertainties. This might indicate that Eddington bias 
(i.e. uncertainties on measurements) could not be sufficient to 
explain the anti-correlation we observe, even if this assumption 
is difficult to assess. The discrepancies between the temperature 
maps derived with a fixed and free pc, with temperature distribution 
not correlated with potential dust heating source in non-barred 
galaxies and the well-known degeneracy between temperature and 
emissivity index of course question the physical sense of letting 
both temperature and emissivity vary at the same time. We also 
caution the fact that even with the good resolution of Herschel 
instruments, the ISM elements studied in a 18'' X 18'' pixel are 
large for the galaxies of the sample, even for our closest object 
(~0.3 kpc for NGC 7793). We thus probe a large mixture of dust 
temperatures for each ISM element. It can therefore be challenging 
to disentangle a real relation or an anti-correlation created by the 
addition of several dust temperatures, with an "artificial" flattening 
of the slope when cold dust is present along the line of sight. We re- 
fer to |Juvela & Ysard| (|2012^ for a detailed analysis of the possible 
eff'ect of temperature mixing on the T - /? relation. A compromise 
to avoid these degeneracies for our non-barred galaxies could be 
to use /5 as a free parameter to derive an average emissivity index 
across the galaxy and then fix p to this value to derive temperature 
and mass surface density maps (as done in [Planck Collaboration] 
[201 lc| for instance). Studies on resolved structures in the Milky 
Way (e.g. |Paradis et al.|2009| ) or in closer galaxies like the Magel- 
lanic Clouds are necessary to restrict the dust temperature range 
probed per resolved elements and contribute to resolve these issues. 



5 DUST MASSES 
5.1 Dust mass maps 

Since most of the dust mass of a galaxy resides in its cold phase, 
we can directly derive dust mass maps from the cold temperature 
maps previously generated. To do so, we use the relation: 

Mdust = — p (8) 

with Sy the flux density, the reference frequency, Ky^ the ab- 
sorption opacity of dust grains at that reference frequency, D the 



distance to the galaxy and B(Tc) the Planck function for a given 
temperature Tc. Dale et al. ( 2012 ) derive global dust masses of the 
KINGFISH sample from MBBs using the reference frequency 
- c I 250 yL/m and the SPIRE 250 //m flux densities. Similar tech- 
niques have been used to derive dust mass surface density maps 
in very nearby galaxies by [Foyle et aT] { [2012 ) for M83, |Parkin| 
|et al.||2012 ) in Centaurus A or Mentuch et al. (submitted to ApJ) 
for the Whirlpool system. Since the 500 jim band is less sensi- 
tive to temperature variations than the 250 yt/m band, we prefer to 
use Vo = c / 500 jim as the reference frequency. We note that us- 
ing 500 jim as a reference wavelength for dust calculation could 
make our masses dependent on submm excess if present in our 
objects. However, we do not expect submm excess to be signif- 
icant at 500 iim for our selection of mostly metal-rich galaxies 
fPraine et ar i |2007l [Gordon et aL||20T0| [Galametz et"ar| [20TT]). 
We choose a dust opacity /c^^X250 yL/m)=3.98 cm^ g~^ (taken from 
www.astro.princeton.edu/~draine/dust/dustmix.html and a Milky 
Way Rv=3.1 model) and for each case {fi fixed or free) derive the 
Kabsi'^^^ yum) such as! 



, 500 



(9) 



In the Milky Way Ry=3.1 model we use, is ~2 in the 100- 
1000 yL/m range. Technically, our choice of fixing the scaling term 
of our dust masses from the reference wavelength 250 fim depend 
on where the intrinsic emissivity deviates from this model. |Li &| 
"DraineJ { [200 1| ) made small modifications to their amorphous sili- 
cate opacity to better match the average high Galactic latitude spec- 
trum measured by COBE-FIRAS from A>250 yt/m. We thus choose 
to modify our emissivity from the same wavelength. We show the 
dust mass surface density map obtained for the galaxy NGC 3627 
with and Pc free (in M© kpc~^, log scale) in Fig.js] We over- 
lay MIPS 24 fim contours for comparison. We gather the dust mass 
surface density maps (with /?c=2.0 and free) for the whole sam- 
ple in appendix in Fig. |A1| (right panels for each galaxy). We ob- 
serve that peaks of this distribution are usually located in the center 
of the galaxies (NGC 337, NGC 3351 for example). Nevertheless, 
peaks are also found at the bar ends in NGC 3627, coinciding with 
peaks in the molecular hydrogen detected through CO observations 
( [Leroy et"aL]|200 9 ), linking the dust concentrations with the cold 
gas reservoirs. This result is similar to that ofl Foyle et al.| ( [20l"2| 
who find that the highest dust concentration in M83 is located in the 
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Table 4. Global dust masses in units of 10^ solar masses 



NGC 337 
NGC 628 
NGC 1097 
NGC 1291 
NGC 1316 
NGC 1512 
NGC 3351 
NGC 3621 
NGC 3627 
NGC 4826 
NGC 7793 



Global fit 
. = 2.0 Pc free [D12] ^ 



361+g 
441^20 

1390!^66 
283+J4 

170!ji 
462^32 

258!J3 
369+}9 

601^29 

76^4 
140+^ 



205^38 
407+66 

^■^ -52 

1380+}90 

151-? 
312^95 

95+22 

2%rl 



892:i» 



225+^3 
291!37 

904-1' 

293+79 

126+13 
164!]^ 
196!g 



294+^0 242+84 



439!^7 
52!^ 



Sum of the dust masses obtained per pixel 
fic = 2.0 ySc = 1.5 ySc^l.O yScfree [A12b] ^ 



347+46 
381+46 

914+45 

^r^_39 
97Q+28 

358^44 
571^24 
75!^ 
153+1^ 



180!}^ 
207+21 
749+69 

127+23 
125+14 
143+20 
160!!4 

191+22 

317!|9 
40+J 
75+1 



94+12 
113^11 
415+38 
73+12 
72!!' 

103+}^ 

21+4 

3rl 



914+40 

371+^^ 
1420+231 
241+44 

-55 

186^30 
195+34 

283+36 
782^80 



48+5 



240+g (306) 
329+1 (355) 
1241J226 (1338) 

194+106 (422) 
129+2} (183) 
188+35 (231) 
238+34 (283) 
333+58 (4Q9) 
475+g (543) 

59+9 (61) 
131^42 (145) 



1 Total dust masses derived by |Dale et al1j2012) ([D12]) using the |Draine & LiH2007) model to fit the global SED. 

2 The dust mass maps are derived by [A12b] using the |Draine & Li||2007) model applied pixel-by -pixel. The dust masses are derived adding the dust masses 
of pixels that are above our signal-to-noise threshold. The total dust masses they derive (whole galaxy) are given in parenthesis. 



nucleus and at the ends of the bar of the galaxy. We observe peaks 
of the dust mass surface density in the northeastern and southwest- 
em arms of NGC 1512, here again at the end of the bar where the 
dust emission peaks from Spitzer/IRAC to Herschel/SFIRE wave- 
lengths. Unfortunately, no CO observations are available for this 
galaxy that would allow comparisons with the molecular gas dis- 
tribution. We observe high concentrations of dust at the northwest 
and southeast sides of the nucleus in NGC 1316, which correspond 
to peaks of dust emission at PACS and SPIRE wavelengths (see 
the resolved blobs in Fig.[TJ. These dusty structures coincide with 
reservoirs of molecular hydrogen detected in CO ( Lanz etal.|2010|. 
[Horellou et al.| ( [200T| have suggested that this molecular gas could 
have an external origin (accretion of material during merging pro- 
cesses). We however remind the reader that NGC 1316 possesses a 
low-luminosity X-ray active galactic nucleus (AGN). Our submil- 
limeter wavelengths thus potentially contain non- thermal emission 
from the AGN that could affect the SEDs and the dust mass map 
we obtain, as discussed in |Bendo et al.| ( |2012f for M81 or in |Parkin| 
|et al.| ( [20T2| ) for Centaurus A. A refinement of our resolved fit com- 
bining non-thermal emission from the AGN and dust associated 
with the X-ray bubbles would be necessary to understand how this 
influences the distribution of the dust mass surface density in this 
object. 

We note that a non-negligible amount of dust is detected in 
the ring of NGC 1291, with higher dust mass surface densities than 
in the nucleus in the west arm in particular. | van Driel et al.|(r988| ) 
observations of Hi in this galaxy have revealed a concentration to- 
ward the ring and a central hole in the distribution. Dust thus seems 
to trace the atomic gas in this galaxy. A proper comparison of the 
ring/nucleus temperature and masses in NGC 1291 is detailed in 
Hinz et al. (submitted to ApJ). We finally note that the choice of fic 
(fixed or free) does not affect the locations of the dust mass surface 
density peaks in our galaxies but directly influences the total dust 
masses derived for each galaxy, as shown in the following section. 



5.2 Dependencies of the derived total dust masses 

In this final section, we aim to quantitatively compare the total dust 
mass derived for our objects with the different modelling methods 
used in this paper. Table [4] summarizes the total dust masses ob- 
tained by adding the individual dust masses determined pixel-by- 
pixel with our two-MBB modelling technique in the four different 
cases j3c=2.0, 1.5, 1.0 or free. The global fit section indicates the 
total dust masses calculated from temperatures derived with a two- 
MBB fit of the integrated fluxes. The errors on the dust mass are 
directly deduced from the temperature errors, the main source of 
uncertainties for our the dust estimates. The non-linearity of the 
Planck function at IR wavelengths is translated into the asymmet- 
ric errors of our dust mass estimates. For the local MBB fits, the 
final errors are obtained by summing the individual uncertainties. 
[A12b] perform a pixel-by-pixel SED modeling of the KINGFISH 
galaxies using the |Draine & Li| ( [2007] ) model and produce dust mass 
maps for the whole sample. We sum up the local dust masses they 
obtain in the regions we selected (3-cr detections in all Herschel 
bands) in our galaxies and add these values to Table[4]for compari- 
son. 

5.2.1 Influence of the emissivity index 

Emissivity spectral indices ranging between 1.0 and 2.0 are com- 
monly used in the literature. Table [4] highlights the strong depen- 
dence of the total dust mass estimate with the choice of pc. Figure|9] 
directly illustrates the decrease of the dust mass with the emissivity 
index. We note for instance that MBB fits using a cold dust emis- 
sivity index of 1.5 lead to a decrease of the total dust mass of ~50 
% compared to those performed with a standard emissivity index of 
2.0 and that MBB fits using an emissivity index as low as 1.0 will 
lead to a decrease of -70 %. For NGC 337 and NGC 7793, the dis- 
tributions of dust temperatures and emissivity indices derived from 
the Monte Carlo simulations on integrated fluxes (see Fig. [3]) seem 
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Figure 9. Dependence of the total dust mass with the cold dust emissiv- 
ity index. The x-axis indicates the emissivity index in the modelling. Dust 
masses are normalized to those deduced with a standard emissivity index of 
2.0. 



to rule out the use of a single temperature dust model with /?c=2.0 
to model the cold dust emission at global scale. At local scales, the 
cold dust emissivity maps of NGC 337 and NGC 7793 show that 
is always lower than 2.0, leading to the same conclusion. We 
note that a MBB fit with /?c=2.0 is not a valid assumption at global 
scale for NGC 1512 either. The emissivity map does not, however, 
show such a trend at local scales. We conclude that the choice in 
the emissivity index value has a direct and significant impact on the 
dust mass deduced and thus on the interpretation made from these 
values (gas-to-dust mass ratios for instance). It also complicates di- 
rect comparisons between studies that call for diff'erent models and 
assumptions. 



5.2.2 Influence of the SPIRE data 

Several studies have stressed the importance of observations long- 
ward of 160 iim to properly determine dust temperature and masses 
( [Gordon et al.|2010| ). |Galametz et aL| ( |2011| show that metal-rich 
galaxies usually have dust SEDs that peak at long wavelengths and 
require constraints above 160 fim to position the peak and sam- 
ple the submm slope of their SEDs and thus correctly probe the 
dust temperature distribution. On the other hand, metal-poor galax- 
ies modeled with submm data could show higher dust masses than 
when modeled without submm data, due to the flattening of the 
submm slope or submm excess often detected in those objects. We 
can directly compare our dust masses derived with submm data and 
those derived with data up to the 160 jjm constraint of Spitzer. We 
observe that dust masses obtained using Herschel data can diff'er 
from those derived with Spitzer-bdi^^d fits available in Draine et aL| 
([2007}, especially for the galaxies NGC 337, NGC 628, NGC 1097 
and NGC 1316. Do we see any trend with metallicity? Most of 
the galaxies of the sample show a decrease in the total dust mass. 
NGC 1316, the most metal-rich galaxy of the sample, shows the 
second highest decrease in our dust masses. The total dust masses 
of NGC 3621 or NGC 7793 - galaxies that have lower metallicities 
- are, on the contrary, rather unaff'ected by the addition of Herschel 
data. Of course, any statistical interpretation of metallicity eff'ects 
is rather hazardous for this small sample. Nevertheless, these trends 
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Figure 10. Dependence of the total dust mass with submm data included in 
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with metallicity are consistent with the dependence expected from 
Galametz et al.|(|2011 ) and are already mentioned in the study of 
Dale et al.| ( |2012| ) performed for the whole KINGFISH sample. 
To test the influence of SPIRE data on the dust mass we obtain 
at global scale, we apply our two-MBB fitting procedure on data 
up to 250 yL/m or data up to 350 fim. The emissivity index is fixed to 
2.0, 1.5 or 1.0. Figure [T0| shows how the total dust mass estimates 
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evolve with the number of SPIRE submm data included in the fits. 
We observe that two MBB fits with /3c=2.0 performed with data 
up to 250 yum could diff'er from those derived using the complete 
SPIRE coverage by +10% (-17% for NGC 1512 due to the very flat 
submm slope of this object not sampled when we restrict the fit to 
data below 250 yum). [A 12a] performed comparisons between dust 
masses derived with data up to 160 yt/m or data including the full 
SPIRE coverage for NGC 628 and NGC 6946 and found similar 
values (diff'erence of ~ 16% and 10% respectively). This diff'erence 
is even more pronounced with ySc=1.5 or 1, reaching up to 30 % in 
the second case. We finally note that using a two-MBB fit and data 
up to 350 yum lead to a systematic underestimation of the dust mass 
(by up to 10%) if fic is fixed to 2 for the galaxies of our sample. 



5.2.3 Comparisons with \Skihha et al.\{2011\ 

ISkibba et alTI ( [20TT] ) have estimated dust masses for KINGFISH 
galaxies using single modified blackbody fits with an emissivity 
index /5 of 1.5 and, as in this study, SPIRE 500 jim flux densities 
as a reference wavelength in their dust mass calculation. From our 
global two-MBB fits with Pc fixed to 1.5, we obtain dust masses 
that are larger by a factor of 1.1 to 2.6 than those obtained with 
the same cold dust emissivity index by |Skibba et al.| ( [2011| ). This 
is probably linked with their use of a unique modified blackbody 
to fit the 70-to-500 yum range that leads to systematically higher 
temperature estimates for the cold dust (warm dust contributing to 
the 70 jim flux as mentioned in §3.4) and thus to lower dust masses. 



5.2.4 Comparisons between MBB fits and the \Draine & Li\ 
( |2007| model 



We observe that, as already noticed in Dal e et al.|(2012| , global 
dust masses derived with the D raine & Lir i2007l model only 
diff'er from models using MBBs and a standard emissivity index 
of 2.0 by less than 35%. A local estimation of the dust masses 
leads to similar results, namely an agreement within 30%. It can 
be explained by the fact that in the wavelength range covered by 
SPIRE bands, the SED of the diff'use ISM component predicted by 
the |Draine & Li| ( |2007 l model (corresponding to the bulk of the 
dust mass in these models) is quite similar to a single temperature 
modified blackbody with J3c = 2. We nevertheless note a significant 
diff'erence between the dust mass obtained with a global fit for 
NGC 1512 when = 2.0 and that obtained by |Dale et al1 ( [20T2l ). 
As mentioned previously, this galaxy shows a flat submm slope 
(global j3c < 1.0, Table|3]). By fixing fie to 2.0, we implicitly let our 
two-MBB model invoke colder dust to account for the SPIRE 500 
jam emission. We observe that the [Draine & Li|(2007|) mo del used 
to fit the global SED of this galaxy in [Dale et al.| ( |20T2] ) tends to 
under-predict the SPIRE 350 and 500 fim fluxes for this galaxy, 
which probably leads to a underestimation of the cold dust mass in 
this particular object. It is thus not surprising that we derive high 
dust masses for this galaxy. This could imply that cold dust could 
be present in this galaxy (and thus more dust mass), especially in 
the outer regions of the galaxy (not modeled in this study if the 
detection is below our chosen signal-to-noise threshold) or that the 
global emissivity index is not well reproduced by standard graphite 
grains used to model carbon dust in the [Draine & Li| ( [2007] ) model. 



5.2.5 Resolution effects 

We want to compare our global dust masses with our resolved dust 
masses. In our local study, we restrict ourselves to regions above 
our signal-to-noise threshold in all Herschel bands in our galaxies. 
Unfortunately, we are consequently missing dust mass in the faint 
outer parts of the galaxies. Making a direct comparison with the 
total dust mass derived from integrated fluxes and the sum of the 
individual dust masses obtained in the diff'erent pixels is thus not 
directly possible from this study. Nevertheless, we observe that 
the dust masses obtained globally using the D raine & Li| ( [2007| ) 
model by [Dale et aL| ( |2012| ) are systematically lower that those 
derived by summing the local dust masses (Table |4] last column, in 
parenthesis), by 11 % for NGC 7793 to 69% for NGC 3621, one 
of the lowest metallicity of the sample. Since dust masses derived 
using a MBB technique and an emissivity index of 2.0 are very 
close to those obtained using the [Draine & Li| ( [2007| ) model, we 
can easily predict the same eff'ect for studies using MBBs, models 
that are intensively used by the scientific community. Those results 
agree with Gall iano et al.[ ( [201 1) who observe similar resolution 
eff'ects while deriving dust masses in a strip of the LMC. The dust 
masses they derive from integrated fluxes are significantly lower 
than dust masses obtained from a local scale study. The dilution 
of cold massive regions in hotter regions is invoked as a possible 
explanation of this eff'ect. If this hypothesis is valid, we should 
expect an increase of our dust mass estimates until the resolution 
of our submm maps reach scales where cold dust dominates (size 
of giant molecular clouds?). Observations at higher resolution, 
with ALMA for instance, will help to address this issue. These 
diff'erences of total dust mass estimates with resolution highlight 
the necessity of a resolved study to properly quantify the total 
dust mass budget of galaxies, when it is possible. Global dust 
mass of high redshift galaxies might be underestimated because 
of the lack of resolution. Moreover, our sample does not cover 
the full diversity of the SEDs of distant galaxies, especially the 
very luminous IR galaxies detected at high-redshift. In some cases, 
resolution eff'ects may be even more significant. The choice of the 
fitting method therefore becomes decisive. 



6 CONCLUSIONS 

We combine SpitzerlMWS and Herschel data to derive global prop- 
erties of the cold dust in a sample of 1 1 nearby galaxies using a 
two-temperature fit to model their SEDs. At global scales, 

• Our objects show a wide range of SEDs, with a clear evolution 
of the lOO-to-24 jjm flux density ratios from actively to quiescent 
star-forming objects and variations on the global cold dust param- 
eter pc we derive, from pc = 0.8 to 2.5. 

• We perform a Monte Carlo simulation in order to quantify the 
efl'ect of uncertainties on the parameters derived from our fitting 
technique. We show that even with an exhaustive sampling of the 
global thermal dust emission, an artificial T-p anti-correlation can 
be created by uncertainties in the flux measurements. 

We use the unprecedented spatial resolution of Herschel to 
derive resolved dust properties within our objects and study the ro- 
bustness of our result with our SED model assumptions (emissivity 
index, SED coverage, resolution). Our study at local scales shows 
that: 

• When we fix the cold emissivity index to a standard Pc-2.0, 
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we obtain a smooth distribution of the cold temperatures, with a 
radial decrease toward the outer parts of galaxies. In some objects, 
the temperature distribution follows the distribution of star-forming 
regions (NGC 3627 for instance). Fixing to 1-5 or 1.0 influences 
the temperature range derived but does not strongly modify the 
temperature relative distribution. 

• When is free, barred spirals (e.g. NGC 1097, NGC 1316, 
NGC 1512, NGC 3351 or NGC 3627) show temperature distribu- 
tions that are similar to those obtained with a fixed Pc and homoge- 
neous emissivity index values across the galaxies. For non-barred 
galaxies, Pc seems to decrease with radius. We nevertheless obtain 
homogeneous cold temperature distributions for those objects. The 
fact that these temperature maps do not seem to relate to any dust 
heating sources questions the use of a free p in those objects. Ap- 
plying a /?-free model to deduce an average emissivity index across 
the galaxy and derive the temperature map using this median value 
could help to avoid potential degeneracies between T and p in those 
objects. 

• Gathering our local results in a Tc-Pc diagram clearly leads to 
a correlation between the two parameters. This correlation could be 
linked with 7 j a real emissivity variation of dust grains with tem- 
perature 2) noise eff'ects and uncertainties on flux measurements 3) 
temperature mixing along the line of sight. 

• The dust mass surface density usually peaks in the center of 
galaxies but also seems to coincide with extra-nuclear peaks of star 
formation at the end of the bars (usually corresponding to molecu- 
lar gas reservoirs). 

• The dust mass estimates seem to be aff'ected by resolution, 
with a systematic increase of total dust masses in the "resolved" 
case compared to dust masses globally determined. 

• The submm coverage of our SEDs with Herschel observations 
modifies the dust masses compared to those derived from Spitzer- 
based fits. Two MBB fits performed with data up to 250 jim diff'er 
from those derived using the complete SPIRE coverage (data up to 
500 iim) by up to 17% if /?,=2.0, up to 30% if /?,=1.0. 

• Dust masses derived globally or on a local basis using the 
[Draine & Li| ( |2007| ) model only diff'er from models using modi- 
fied blackbodies and Pc=2.0 by less than 30%. Dust masses derived 
from modified blackbody fits using a shallower parameter p^ are 
systematically lower that those derived with Pc=2.0 (lower by 25 to 
46 % if p^l.5 for instance in our sample). Our results highlight the 
necessity to better investigate the variations of the cold dust proper- 
ties to properly quantify the total dust mass if we aim for instance, 
to use it as a gas tracer. 



This study aimed to investigate the properties of cold dust un- 
veiled by Herschel observations of nearby galaxies. Questions re- 
main due to possible degeneracies between the grain emissivities 
and temperatures. Moreover, this paper uses isothermal fits, which 
implies that the intrinsic emissivity index of the dust grains and 
the emissivity index directly obtained from our fits are one and the 
same. The presence of several dust grain temperatures in our ob- 
jects both locally and at global scales also makes a direct extrap- 
olation of dust properties difficult. Nevertheless, combining those 
results with other derived quantities such as the elemental abun- 
dances in the ISM or reliable gas-to-dust mass ratios as well as 
further investigations on closer objects will help us to disentangle 
between the different hypotheses. 
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APPENDIX A: DUST TEMPERATURE, EMISSIVITY 
INDEX AND DUST MASS SURFACE DENSITY MAPS 
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Figure Al. For each galaxy, North is up, East is left. Top: MIPS 24 //m image (MJy sr~^), cold dust temperature map (in K) and dust mass surface density 
map (in M© kpc~^, log scale) derived with /?c = 2.0. All images are convolved to the SPIRE500 PSF before deriving the parameter maps. The FWHM of the 
PSF of SPIRE 500 //m is shown by the dashed circle in the top middle panel for each galaxy. Bottom: Emissivity index map, cold dust temperature map (in K) 
and dust mass surface density map (in M© kpc~^, log scale) derived with T^. and Pc as free parameters. The value ySc=2.0 is indicated on the color scale of the 
emissivity index map. We overlay the MIPS 24 //m contours (smoothed to the resolution of the temperature and beta maps) for comparison. 
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Figure Al. Continued. 
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Figure Al. Continued. 
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Figure Al. Continued. 
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Figure Al. Continued. 
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Figure Al. Continued. 



